Pair Density Wave Instability of Cold Fermionic Atoms in an Optical Lattice 
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We show that the band insulating state of cold fermionic atoms in a three-dimensional cubic optical 
lattice is unstable to forming a pair-density wave (PDW) state in the presence of strong attractive 
interactions induced by proximity to a Feshbach resonance. The PDW state in this model is a 
special type of supersolid, which arises from the condensation of predominantly interband Cooper 
pairs at a nonzero, generically incommensurate, wavevector, and is expected to be accompanied by 
an induced density modulation. We also discuss the mean field phase diagram of a two-band tight 
binding model, which exhibits a PDW state, and discuss the relation of the PDW instability in this 
model to the Halperin-Rice exciton condensation instability in indirect bandgap semiconductors. 



Introduction. — The theme of coexisting or com- 
peting order parameters is common to several strongly 
correlated systems, including high temperature cuprate 
[l[ and pnictide Q superconductors. Most notably, there 
have been striking experimental indications of stripes or 
checkerboard patterns of spin and charge modulations, 
that coexist with superconductivity in several cuprate 
materials 0, 0, 0, @, 0j- Motivated by the observa- 
tion 0] of quasi-two-dimensional superconductivity co- 
existing with stripe order in the layered superconductor 
Lai.87sBao.i25Cu04 (LBCO), Berg et al. [9|] have pro- 
posed, on phenomenological grounds, that a distinct state 
of matter, named a 'pair density wave' (PDW), is realized 
in this material. In its simplest avatar, the PDW state re- 
sults from condensing singlet Cooper pairs with nonzero 
center-of-mass momenta ±Q and is accompanied by an 
induced charge density modulation at momenta ±2Q. It 
is thus similar to the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state of magnetized superfiuids [lfj except that 
the PDW does not require broken time-reversal symme- 
try. In contrast to earlier proposals for a Cooper pair in- 
sulator 111, [l2j in LBCO, the PDW state is a supersolid, 
in that it has coexisting superfluid and density orders, 
which break lattice symmetries. However, it is very dif- 
ferent from the supersolid state proposed to exist in 4 He 
[l3L [p| , or the supersolids realized in simple lattice bo- 
son models [HI, since the superfluid order parameter in 
the PDW state has no uniform Fourier component. The 
bosonic analog of the PDW may be found in lattice mod- 
els of bosons in which the kinetic energy is 'frustrated' so 
that the single-particle dispersion has multiple minima, 
at nonzero momenta, into which bosons can condense (lfl | 
without breaking time-reversal symmetry. 

In this Letter, we show that there are simple micro- 
scopic models of fermions with short-range attractive 
interactions which support PDW ground states. Our 
work complements earlier Landau theory descriptions 
and Josephson junction models of the PDW state 
Our first example is a one-channel model of fermionic 
atoms near a Feshbach resonance 17, Il8f confined to 



a cubi c op tical lattice. It has been demonstrated re- 
cently [19( that this system shows a superfluid to band- 
insulator transition 20, 21 , 2^ when the lattice depth is 
varied at a commensurate density of two atoms per lat- 
tice site. Here we show, via a more careful study, that a 
PDW state is expected to intervene between the uniform 
superfluid and the band insulator. Our second exam- 
ple is a two-band tight binding model where the chem- 
ical potential is chosen to be such that both bands are 
empty. An appropriate choice of local attractive inter- 
actions between the fermions is then shown to lead to a 
PDW instability of this zero-filling band insulator. 

The key physics which leads to the emergence of the 
PDW state in both these models is the presence of multi- 
ple bands and the dominance of interband Cooper pair- 
ing. In the cold atom model, we present arguments to 
show that, in contrast to intraband pairing, the phase 
space for interband pairing is expanded at nonzero pair- 
ing momenta, which stabilizes an incommensurate PDW 
state. In the two-band tight binding model, the reason 
for the occurrence of the PDW state is that the lowest en- 
ergy momentum points in each band differ by a nonzero 
wavevector Q, which leads to a large Cooper pair sus- 
ceptibility at this wavevector. For this model, we present 
the mean field phase diagram and show that the PDW 
instability is closely related to the Halperin-Rice exciton 
condensation instability in indirect bandgap semiconduc- 
tors [i^l , and some models of Iron-pnictides [24[ . 

Cold Atoms in an Optical Lattice. — We describe 
fermionic atoms with attractive interactions in a periodic 
potential 2^, 21, 2^ using the Hamiltonian (H = 1): 



H = 



d 3 r 



2m 



T c l c i c T 



(1) 



Owing to universality in the unitarity regime, this simple 
theory provides a faithful description of fermionic cold 
atoms tuned near a broad Feshbach resonance [25|. We 
will study this model using mean-field theory which is 
known to be a reasonable approximation near unitarity 
for the qualitative points we wish to make. Fluctuations 
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can be treated systematically using, for example, large- N 
2ll 25L [2^ | , but we will not pursue this here. 



expansions 



We work with a simple cubic lattice potential: 



V r = V 
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where or, is lattice spacing. The quantum numbers 
of single-particle Bloch eigenstates in this potential are 
crystal wavevector k = (k x , k y , k z ) inside the first Bril- 
louin zone (BZ) —ir/aL < k x ,k y ,k z < n/a^, and band 
index n = (n x , n y , \ 



We label the Bloch wavefunctions 



by V'n.fc( r ) an( i the corresponding energies by e n k- 

Near unitarity, the cutoff-dependent contact interac- 
tion parameter U is related to the scattering length a: 
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(3) 



Band- index cutoff, discussed below, is implicit in ([3]). 

A T = superfluid-insulator transition for an even 
number of fermions per site occurs at a critical value of 
the lattice amplitude V, which is a universal function of 

I a and the fermion density 2^, 21 , 2^ • Starting from a 
band-insulating state, the onset of pairing in the mean- 
field approximation can be extracted from the inverse 
static pairing susceptibility matrix: 



n, 



Gq.G'q' 



E 



x r 



n 1 k 1 ;n 2 k 2 n\k\\n 2 k 2 



d 3 fci d 3 fc 2 /(Ufcx)-/(-Ufc a ) 
(2tt) 3 (2tt) 3 

i2z)3 -6(q-q')6 GG >, (4) 
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where q are first BZ wavevectors, G are reciprocal lattice 
vectors, = e n k — /i, /(C) is Fermi-Dirac distribution 
function and Y are vertex functions: 



rc-ifci :n 2 k 2 



= / d J r$^(r)^ ifci (r)^ 2fc2 (r). (5) 



Since crystal momentum is conserved, = 
IIgg'(9) x (2tt) 3 ^(<Z — q'). We will use the plane wave 
representation for TlcG'(q), corresponding to the pair 
wavefunctions $g<j(t ) = e^ q+G ' r . All eigenvalues of the 
matrix T1gg> are positive in the band insulating state. 
When the lowest eigenvalue n(q) becomes negative at 
some wavevector q = Q, the insulating state becomes 
unstable to a superfluid of fermion pairs condensing at 
momentum Q, which is a PDW state 0] if Q 7^ 0. 

Normally one would expect the pair condensation to 
occur at Q = 0. This is certainly true in any single-band 
model of lattice fermions. However, as we demonstrate 
below, interband pairing in multi-band models can give 
rise to pairing instability at a finite Q. Figure Q] shows a 
contour plot of 11(g) for two fermions per lattice site at 
T = 0. Each contour in this plot j2?| defines a superfluid- 
insulator transition for a different scattering length a. 
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FIG. 1: Contour plot of the lowest eigenvalue of the inverse 
pairing susceptibility 11(g) for q — (q,q,q), indicating pairing 
instability for different values of scattering length a. PDW 
ordering wavevector Q at the onset of transition is plotted 
with a bright solid line as a function of inverse lattice depth 
(E r = it 1 1 Amah is molecular recoil energy). 



Insulators are obtained in the limit of large V, so the 
largest V on the contour for a given a is the critical lat- 
tice depth V c (a) for the onset of pairing. One can see 
that V c (a) generally corresponds to a nonzero value of 
q = Q. The smooth evolution of \Q\ with the lattice 
depth in the deep BCS limit indicates that the formed 
PDW state is incommensurate. As the pairing interac- 
tions become stronger moving toward the BEC limit, \Q\ 
grows and possibly eventually saturates at the BZ edge 
making the PDW commensurate. Another notable fea- 
ture is that sometimes a superfluid at not too small q 
can be destabilized by both increasing and decreasing V 
(for example, the dashed contour in FigfT]). The latter 
illustrates that interband pairing is responsible for su- 
perfluidity at finite q, which can be expected to weaken 
with decreasing V. Note that Mott physics plays no role 
in this instability, since repulsive interactions between 
Cooper pairs are captured only by two-boson correlation 
functions. Also, without knowing such correlations we 
cannot rule out time-reversal symmetry breaking. 

Even though the PDW instability we observe is a weak 
feature, we found its robust convergence with increasing 
cutoffs. In the following we provide simple arguments to 
show why multiband models can favor a PDW instability. 

The incommensurate PDW owes its existence to phase- 
space restrictions for interband pairing. Consider two 
particular fermion bands along some momentum direc- 
tion in the first BZ, separated by an indirect 'gap' (which 
may be filled by other bands). For concreteness, let us de- 
scribe these two bands using a one-dimensional toy model 
with cil = 1 illustrated in Fig. O Since momentum is 
conserved only modulo reciprocal lattice vectors G, we 
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FIG. 2: Pairing of two fermions with crystal momenta ki and 
&2- Intraband pairing in (a) occurs when ki + k2 — q = 0. 
Interband pairing in (b) occurs when ki + k^ — q = — 2tv 
(assuming q > 0). Thick arrows show the trajectories of fci 
and k<2 dictated by momentum conservation. 



can rewrite the vertex functions 



as: 



X) < 9 ™ 2 (G') x 27r<5(fc 1 + k 2 - q + G') , (6) 



where the coefficients (£?') depend on details of the 

band-structure. In the limit of small lattice depth V all 
but one of these coefficients for any fixed (ni, n 2 , G) be- 
come small, and must gradually vanish as V — > 0, when 
momentum conservation becomes exact. The remain- 
ing large coefficient occurs at a value of G' , which de- 
pends on (ni,n,2,G), and approaches unity as V — > 0. 
For example, pairing into a plane- wave superfluid at 
q £ first BZ is given by Al 1 (G) rs 8g,o (intraband) and 
Ai2 ((*) — A\x (G) ~ ^G.27rsgn(g) (interband) for small V. 
We simplify the following discussion by focusing only on 
this pairing channel which reduces the inverse pairing 
susceptibility matrix Hqc (q) to a scalar H(q). The other 
allowed pairing channels, involving condensate harmon- 
ics at larger reciprocal lattice vectors, will be neglected 
since they have relatively small amplitudes at small V, 
and only slightly enhance the PDW instability. 

Using ([6|) we find that the main contribution to intra- 
band pairing for small V at T = comes from 



n (bD( g )~_i. J dkldk , 



S(ki + k 2 — q) 
£i,/ci + £i,fc 2 



which is illustrated in Fig. [2ja). Since k\ and k 2 are 
restricted to the first BZ, the number of states available 
for intraband pairing decreases with q. Consequently, the 
magnitude of IB 1,1 )^) decreases with q and thus purely 
intraband pairing would occur at q = 0. The dominant 
interband contribution 



n 



(1,2) 
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2^ 
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6(ki + k 2 + 27rsgn(q) - q) 
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illustrated in Fig. E^b) has the opposite behavior be- 
cause the number of states available for interband pairing 
increases with q. Therefore, interband processes alone 



would prefer pairs to condense at a BZ edge, i.e. at the 
largest possible value of q. 

It is important to note that H(q) ~ \q\ for q — > due to 
the boundaries of momentum integrals in all n^™ 1 ™ 2 ^ (q), 
as can be seen from Fig. Only for V — and in the 
tight-binding limit do these linear contributions cancel 
out, leading to U(q) ~ q 2 . The initially negative slope of 
H(q) leads to a local minimum at q ^ 0. The location of 
this minimum is determined by the relative strengths of 
interband and intraband contributions, so that in prin- 
ciple it can be anywhere in the BZ, making the PDW 
incommensurate. However, for large lattice amplitudes 
V a tight-binding limit provides a better description of 
the system and naturally yields a commensurate PDW 
as we illustrate below. 

A linear H(q) for q — > is incompatible with a uniform 
superfluid instability. Since phase-space restrictions for 
pairing in the presence of a periodic potential generally 
result in a linear H(q), we argue that a PDW supersolid 
always preempts an ordinary superfluid instability of the 
band insulator. This is consistent with our numerical 
findings. Note that fluctuations beyond the mean-field 
approximation cannot remove the PDW instability unless 
some accidental cancellations of the linear terms occur. 

Two-band tight binding model. — Let us next 
turn to a tight-binding model which is of interest for 
fermions in deep optical lattices or for strongly correlated 
solid state materials. We assume a model Hamiltonian 

(i.j)ncr ina 

~ U 53 (^™in2 C L 1 T C m 2 i)(^ 1 £ 2 C 4 f 2 T C ifii)' ( 7 ) 
i,ni ,n 2 ,t\ ,£2 

where the fermions c have a band-index n in addition to 
spin a. The noninteracting ground state of this model is 
either a normal metal or a band insulator, depending on 
the parameters in the Hamiltonian. 

Upon including interactions, we find that the inverse 
Cooper pair susceptibility of the noninteracting ground 
state is given by 



r%) 



1 



U N. 



(f(£n,k)-f(-tje,-k+q)) 



(8) 



where = — 2£„(cos k x + cos k y ) +7„ — and N s is the 
number of lattice sites. In contrast to the previous model, 
the momentum dependence of the susceptibility in this 
tight binding limit arises only from the band dispersion, 
while the coupling constants are momentum independent 
unlike the pairing vertices studied earlier. For simplicity, 
we focus here on a two-dimensional two-band model and 
restrict ourselves to the case where An = X 22 = cos 6, 
and A12 = -A 2 i = sin^, with < 6 < vr/2. With this 
parametrization, the overall pairing strength is controlled 
by U, while tuning the angle 9 takes us from pure intra- 
band pairing (9 — 0) to pure interband pairing (9 = 7r/2). 
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FIG. 3: Mean field phase diagram of the two-band tight- 
binding Hamiltonian discussed in the text with dispersion pa- 
rameters ti = l,ta — — l)7l — 0,72 = 4, and fi — —4.5. Thin 
(thick) lines indicate second (first) order transitions. 



Fig.[3]shows the mean field phase diagram of this model 
for a specific choice of dispersion and chemical potential 
at which the noninteracting state is a zero-filling band 
insulator. We find that this band insulator can undergo 
continuous transitions into either a uniform superfluid or 
a PDW state depending on whether intraband or inter- 
band interactions dominate. The PDW instability in this 
model arises from the fact that the dispersion minima of 
the two bands differ in momentum by Q = (tt, 7r), so 
that there is a divergence of the interband Cooper pair 
susceptibility at this wavevector for sufficiently strong in- 
terband interactions. Upon making a particle-hole trans- 
formation (followed by a spin rotation) of the fermions 
in the lower band (band-1), it is easy to see that the 
interband singlet Cooper pair condensate maps onto an 
exciton condensate, where the exciton is made of a hole 
in band-1 and a particle in band-2, whose momenta dif- 
fer by Q = (w,tt). This PDW instability can thus be 
recognized as the particle-particle analog of the Halperin- 
Rice exciton condensation instability in indirect bandgap 
semiconductors 23]. The superfluid to PDW transition 
is first-order for this model. In other models, or for other 
parameter values, there could be an intervening conven- 
tional supersolid state [281 ]. 

Discussion. — We have shown that the PDW state 
arises in simple and experimentally relevant microscopic 
lattice fermion models. In the multiband tight binding 
limit, the PDW state appears to be robust provided one 
has strong interband interactions. The PDW instability 
for cold atoms near unitarity in a cubic optical lattice is 
a comparatively weaker effect at the lattice depths and 
densities we could explore. Also, the ordering wavevec- 
tor is incommensurate and small, which could make it 
difficult to observe this effect in cold atom experiments. 
However, we expect it to be possible to significantly en- 
hance the PDW modulation by engineering different lat- 
tice potentials and using higher atom densities. 
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